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SECTION  I 


INTRODUCTION 

The  work  reported  in  this  study  concerns  indicia!  unsteady  aerodynamic 
theory  and  its  role  in  conventional  aircraft  dynamics.  The  basic  objective 
of  this  effort  was  to  investigate  the  extent  to  which  improved  knowledge  of 

the  unsteady  derivatives  C^  ,  Cm  ,  etc.  would  allow  for  extension  of  the 

a  a 

boundaries  which  constrain  the  full  potential  operation  of  existing  high 
performance  aircraft.  To  accomplish  that  task,  it  was  found  necessary  to 
develop  a  simplified  but  accurate  theory  for  the  indicial  unsteady  aerodynamics. 
The  simplified  concepts  developed  in  this  report  have  been  extended  to  a  wider 
range  of  wing  geometries  and  considered  for  application  to  the  parameter 
estimation  problem  as  reported  in  references  [1]  -  [2].  The  results  of  the 
analysis  are  presented  following  a  brief  overview  of  classical  and  contemporary 
unsteady  aerodynamic  theory. 


SECTION  II 

REVIEW  OF  UNSTEADY  WING  THEORY 
1 .  Classical  Research 

A  systematic  study  of  unsteady  aerodynamics  began  shortly  after  Prandtl's 
theory  of  bound  vortices  was  published  in  1918.  This  was  followed  by  an  analy¬ 
sis  of  Birnbaum  [3]  in  which  he  calculated  the  aerodynamic  force  on  a  harmon¬ 
ically  oscillating  two  dimensional  flat  plate.  One  of  the  more  referenced  and 
fundamental  studies  was  made  by  Wagner  [4]  in  which  he  calculated  the  unsteady 
lift  on  an  airfoil  impulsively  started  from  rest  or  equivalently  an  abrupt 
change  in  angle  of  attack  of  a  plate  in  a  uniform  stream.  Glauert  [5]  extended 
Wagner's  method  to  the  problem  of  computing  the  forces  and  moments  on  an  oscil¬ 
lating  airfoil.  In  1935,  Theodorsen  [6]  published  his  exact  solution  for  the 
lift  on  a  harmonically  oscillating  airfoil  with  a  flap.  This  work  introduced 
the  Theodorsen  F  and  G  functions  to  unsteady  aerodynamics  theory.  Kiissner  [7], 
more  or  less  closed  this  creative  era  with  his  solution  for  the  lift  on  an  air¬ 
foil  as  it  penetrates  a  sharp  edged  gust.  The  mathematical  functions  associated 
with  the  works  of  Wagner,  Theodorsen,  and  Kiissner  were  shown  to  be  fundamentally 
related  by  Garrick  [8]. 


During  the  next  two  decades  these  fundamental  solutions  were  extended 
and  sharpened  by  another  generation  of  researchers  [9-16].  Von  K£rm£n  and 
Sears  [9]  developed  a  uniform  circulation  theory  which  treats  the  previous 
studies  as  special  cases.  The  KSrmSn-Sears  theory  is  outlined  in  more  detail 
in  the  Appendix.  Jones  [10]  used  the  basic  results  of  Wagner  to  solve  for  the 
unsteady  lift  on  a  wing  of  finite  aspect  ratio.  Another  significant  contribu¬ 
tion  in  Jones'  work  was  the  development  of  an  expression  for  the  start-up  lift 
on  the  wing.  Further  work  was  performed  on  the  finite  wing  by  Reissner  [12] 
in  which  he  formulated,  in  terms  of  an  integral  equation,  and  solved  for  the 
lift  distribution  on  oscillating  wings.  Extension  of  these  results  to  include 
compressibility  effects  were  performed  by  Lomax,  et  al.  [15]  for  two  and  three 
dimensional  wings  at  high  speeds. 

2.  Contemporary  Research 

The  direction  of  unsteady  lifting  surface  theory  during  the  past  decade 
has  been  toward  digital  computer  application  to  the  classical  formulations 
[17-25].  Ashley,  et  al .  [17]  started  this  movement  in  a  research  article  in 
which  he  proposed  the  use  of  modern  digital  computers  to  obtain  an  exact  solution 
to  the  oscillating  wing  problem.  In  particular,  he  outlined  the  method  now 
referred  to  as  the  "doublet  lattice  method"  for  numerically  solving  the  singu¬ 
lar  integral  equation  for  the  normal  velocity  at  the  surface.  Albano  and 
Rodden  [19]  developed  the  method  to  an  operational  form  adaptable  to  complex 
wing  geometries.  The  method  is  outlined  in  detail  in  Reference  14.  The  basic 
idea  is  to  divide  the  wing  surface  into  trapezoidal  panels.  The  1/4  chord  line 
of  each  trapezoid  contains  a  distribution  of  acceleration  potential  doublets  of 
uniform  but  unknown  strength.  The  strength  of  the  doublets  is  then  determined 
by  use  of  the  flow  tangency  condition  at  the  control  points  of  each  trapezoid 
of  the  wing. 

A  vortex-lattice  method  of  determining  the  transient  lift  buildup  on  a 
finite  wing  has  been  presented  by  Belotserkovskii  [21].  This  method  has  been 
further  developed  by  Atta  [22]  to  account  for  the  wake  geometry  as  well  as  the 
wing-tip  vortex  system.  Due  to  the  likely  importance  of  this  particular  numer¬ 
ical  method  in  computing  indicial  aerodynamic  effects,  a  brief  description  of 
the  method  will  be  given  using  Figure  1  which  indicates  the  manner  in  which  the 
lattice  is  formed  in  the  wake  and  on  the  wing  surface.  At  time  t=l,  the  initial 
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free  vortex  with  strength  6  ^  is  shed  in  the  wake  and  is  pljced  at  an  appro¬ 
priate  location  forcing  satisfaction  of  the  Kutta  condition  at  the  trailing 
edge.  At  this  time,  the  strengths  of  the  wing  bound  vortex  are  solved 

for  while  satisfying  the  flow  tangency  co^jtions  at  the  control  points.  At 
time  t  =  2,  the  free  vortex  of  strength  6  .  is  carried  downstream  while  its 

strength  remains  constant.  Another  free  vortex  with  strength  6^  forms  where 

(1 )  ( 1  (21 
6  ^  was.  Again,  the  strengths  of  the  entire  vortex  field  nv  and  6'.'  are 

solved  for  from  the  flow  tangency  conditions.  This  process  is  repeated  until 

no  further  change  occurs  in  r(^  at  which  time  steady  state  conditions  are 


reached. 

Morino  [23]  has  developed  a  numerical  procedure  to  solve  for  the  unsteady 
lift  on  a  wing  based  on  a  Green's  function  formulation.  The  method  is  appli¬ 
cable  to  either  oscillatory  or  transient  flow  phenomena  and  applies  to  both 
subsonic  and  supersonic  flow.  Basically,  Morino  develops  an  integral  equation 
for  the  unsteady  velocity  potential  which  involves  values  of  the  potential  and 
its  normal  derivative  on  the  surface  of  the  body  and  wake.  He  then  applies  a 
finite  element  method  for  solving  for  the  potential  function. 

The  most  recent  contribution  to  the  field  of  unsteady  aerodynamics  is  due 
to  Edwards  [25].  He  addresses  the  problem  of  active  control  of  supercritical 
flutter  modes.  His  approach  is  to  solve  the  basic  linear  governing  equation 
for  the  velocity  potential  by  use  of  Laplace  transform  techniques. 


SECTION  III 

APPROXIMATE  AERODYNAMICS  ANALYSES 

For  use  in  preliminary  design  and  stability  and  control  analyses,  simpli¬ 
fied  analytical  solutions  are  preferred  to  more  exact  numerical  ones.  Accord¬ 
ingly,  an  approximate  solution  to  the  basic  Wagner  problem  was  sought  which 
preserves  much  of  the  accuracy  of  the  exact  solution  and  physical  principles 
involved.  The  simplified  aerodynamic  model  obtained  is  based  on  the  physical 
principle  of  Prandtl's  lifting  line  theory  and  trailing  vortex  concept.  The 
wake  is  assumed  to  be  compressed  to  a  single  shed  vortex  element  of  appropriate 
strength  moving  downstream  at  a  speed  sufficient  to  approximate  the  Wagner  func¬ 
tion. 
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1. 


Two-Dimensional  Wings 

In  the  simplified  analysis,  the  wing  is  represented  by  a  bound  vortex 
located  at  the  wing  quarter-chord  line  (Figure  2).  The  boundary  condition 
of  no  flow  through  the  wing  is  satisfied  at  the  3/4  chord  line.  A  step 
change  in  wing  angle  of  attack  increases  the  strength  of  the  bound  vortex 
by  an  amount  F.  An  equal  amount  of  vorticity  of  opposite  sign  is  shed 
behind  the  wing  to  form  a  continuous  vortex  sheet  which  increases  in  length 
with  time.  For  the  approximation,  the  shed  vortex  sheet  was  replaced  by  a 
single  shed  vortex  having  the  same  strength  as  the  bound  one  and  moving  with 
a  velocity  KU  as  shown  in  Figure  2.  The  distance  from  the  3/4  chord  point 
at  which  the  vortex  is  shed,  xQ,  is  determined  by  forcing  the  start-up  lift 
to  match  the  results  of  Jones  [10].  For  two  dimensional  wings,  xq  =  c/2. 

The  downwash  at  the  3/4  chord  point  is 


w 


_F 

2tt 


1 

c/2  +  KUt 


] 


(1) 


For  a  unit  step  change  in  angle  of  attack,  w  =  1)  resulting  in 
„  _  ttIIc 


1  + 
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(2) 


2KUt 


Use  of  this  value  of  r  in  the  Kutta-Joukowski  equation  results  in  the  lift 
coefficient 


Ac£  =  2tt 


_} _ 

242KUt 


(3) 


It  was  found  by  comparison  with  exact  results  that  a  value  of  K  =  1/2  resulted 
in  excellent  agreement.  It  is  interesting  to  note  that  this  value  for  Ac^  agrees 
with  an  empirical  value  obtained  by  a  curve  fit  to  Wagner's  solution  by  Garrick 
[8]. 

To  get  the  change  in  lift  coeffecient  for  an  arbitrary  change  in  angle  of 
attack  a(t),  Duhamel's  integral  formula  is  used;  that  is, 


ci(t> 


1 


1 

2  +  U(t-x) 
c  J 


&(t)dx 


(4) 
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To  illustrate  the  use  of  equation  (4) 
a  =  k‘t.  From  equation  (4),  we  get 

cjt)  =  2TTk't~?J^  in  (1  +  Jji) 
To  get  Cj,  differentiate  with  respect 


to  compute  c^.  consider  the  case  where 

a 


(5) 


^  irt  -  4  -rr  in  (1  + 


(6) 


The  hysteresis  effect  due  to  the  unsteadiness  as  computed  from  equation  (5) 
is  shown  in  Figure  3  for  U  =  100  fps  and  c  =  5  ft.  The  curve  labeled  a  >  0 
corresponds  to  a  increasing  linearly  up  to  10°  (k1  =  it/180).  The  angle  of 
attack  is  then  held  constant  until  steady  state  is  reached.  The  curve  labeled 
a  <  0  corresponds  to  a  decreasing  linearly  from  10°  to  0°  (k1  =  -  7r/180).  The 
curve  labeled  "quasi-steady"  corresponds  to  the  instantaneous  lift  buildup  due 
to  a  finite  change  in  angle-of-attack. 


2.  Three-Dimensional  Untapered  Wings 

The  Wagner  function  for  3-dimensional  wings  can  be  approximatea  utilizing 
the  same  concepts  as  for  2-dimensional  wings  as  shown  in  Figure  4  for  a  straiqht 
wing.  At  the  control  point  P,  a  unit  change  in  a  results  in  a  downwash  of 
val  ue 


Aw  =  ^  [  ’V+l  +  ^  +  (2-|)2] 


(7) 


where  b  is  the  wing  span,  A  the  aspect  ratio  and 


2  =  x  +ut 
L  xo  IT 


(8) 


The  corresponding  values  of  circulation  and  lift  coefficient  are  computed  as 


_ rrbU _ 

A2+l  +  /l+(b/2z)2 

2ttA _ 

*AZ+1  +  ^l+(b/2z)2 


(9) 

(10) 


5 


* 


For  an  arbitrary  change  in  angle  of  attack,  the  unsteady  lift  coefficient 
is  computed  by  Duhamel’s  integral  formula  as 


=  2ttA  j 


a(j)  dx 


0  ^7  +  ^ 


(id 


To  find  the  location  of  the  starting  point  xQ,  we  use  the  initial  value  due 
to  Jones  [10]  for  an  elliptical  planform  wing,  i.e.. 


AC^(o)  =  tt/E 


(12) 


where  E  =  Semi-perimeter  =  A+l 
span  A 

This  results  in 


xo/c  = 


A/2 


[(2 (A+l )  -  /aZ+1)Z  -1]I/Z 


(13) 


A  plot  of  x  / c  versus  aspect  ratio  is  shown  in  Figure  5. 


The  lift  as  computed  from  equation  (11)  is  well  suited  for  computation  of 

the  unsteady  stability  derivatives  CL  and  Cm  .  For  instance,  CL  due  to  a 

a  a  a 

linear  change  in  a  is  computed  from  equation  (11)  as  (letting  xQ=  c/2  which  is 
the  limiting  value  for  large  aspect  ratio  wings) 
yTZ - 


r  =  M  A  +  1 

V  c  A 


U 


1  + 


r 

J  t  C 

’A2  +  2  ("I  +  Jp-J- 

T*  i 

L  i/a*  +  l 

*  f 

2^  +  1 

—  All 

2c 


VA2  +  1  +  ^ 


*2 1,  rp 


(14) 


The  Wagner  function  for  the  three-dimensional  swept  wing  can  be  calculated 
by  consideration  of  the  geometry  shown  in  Figure  6.  The  calculation  is  carried 
out  for  a  swept  wing  with  zero  taper  undergoing  a  unit  step  change  in  angle  of 
attack. 
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The  change  in  downwash  at  point  P  is  computed  as 


w  =  [  ^sec2  A  -2A  tan  A  +1  +  ^  ^b2+4z2 


05) 


(16) 


2zb 

+  4  tan  A  (l-§r)  h-  bJinJL 

W 

The  corresponding  change  in  lift  is 

ACl  =  2ttA  [^sec7  A  -2A  tan  A+l  +  ^ 

+  4  tan  A(l-i-)  +  b-  tan  A  ]1/2 

2  2z  /I - 7 

V+4z^ 


3.  Horizontal  Tail  Lift 

In  order  to  calculate  the  unsteady  tail  lift  load,  an  expression  must 
first  be  obtained  for  the  unsteady  downwash  at  the  tail  surface  to  the  wing. 

For  this  calculation,  the  shed  vortex  is  allowed  to  move  downstream  at  the 
free  stream  velocity.  This  assumption  causes  the  initial  disturbance  to  reach 
the  tail  surface  at  the  proper  time.  Numerical  tests  indicate  that  the  time 
variation  in  downwash  at  this  station  is  not  substantially  different  for  either 
assumption  of  shed  vortex  velocity  equal  to  U  or  U/2.  The  geometry  for  the 
calculation  of  downwash  at  the  quarter  chord  point  of  the  horizontal  tail  sur¬ 
face  (point  P)  is  presented  in  Figure  7.  The  downwash  angle  at  point  P  due  to 
the  effect  of  the  wing  wake  is  (for  unit  step  in  wing  angle  of  attack) 


Ac  = 


4rfl  «7  *  /bt4(c/2+X0-VUt)2' 

^  L*t  t  C/2+X  -JL+Ut 


(17) 


where 


r  = 


itUb 


n+A 


+  ^\+ 


f~"b  7 
2(xQ+Ut) 


(18) 


The  effective  angle  of  attack  at  the  tail  surface  (for  small  rotation)  is 
equal  to  the  wing  angle  of  attack  less  the  downwash  angle  a-e.  This  provides 
an  estimate  of  the  tail  lift  for  a  unit  step  change  in  effective  angle  of 
attack  as 


(19) 


where  the  subscript  "t"  refers  to  quantities  based  on  horizontal  tail  surface 
geometry.  For  arbitrary  changes  in  angle  attack,  the  unsteady  tail  surface 
lift  coefficient  is 


CL  (t)  =  AC.  (t-T)[oi(T)-e(T)]  dt  (20) 

t  0  t 


where 


t 

e(t)  =  /  Ae(t-x)  a(x)  dr 
o 


SECTION  IV 

AIRCRAFT  DYNAMICS:  LONGITUDINAL  MODES 


1 •  Equations  of  Motion 

For  a  dynamic  model  of  the  longitudinal  motion,  we  assume  the  following 
perturbation  equations  of  motion 


(21) 

‘■-♦•S', 

(22) 

Cm 

(23) 

6  =  q 

(24) 

The  axial  force  coefficient  for  the  airplane  is  assumed  to  be  a  function  of 
angle  of  attack  only  according  to 
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(25) 


The  normal  force  coefficient  for  the  airplane  is  taken  as  the  sum  of  the 
unsteady  wing  lift,  the  unsteady  tail  lift  and  elevator  contribution  according 
to 


C  —  /  AC.  (t-x )a(x)dx  -  /  AC.  (t-x ) [d(x)-e(x) ]dx  +  C  6  (26) 

o  Lw  o  Lt  z6„  e 

e 

The  moment  coefficient  for  the  airplane  is  taken  as  the  sum  of  the  effects  of 
the  tail  lift,  the  contribution  due  to  pitch,  and  that  due  to  elevator  accord¬ 
ing  to 

*  s  t 

cm  =  -  /0  ACL  (t-x)[d(i)  -  e(x)]dx  (27) 


+  C. 


+  c 

2U  mc 


Reference  to  equations  (21 ) - ( 24 )  indicates  that  u  and  6  can  be  computed  once 
q  and  a  are  determined.  Integro-differential  equations  for  a  and  q  are 
obtained  upon  substitution  from  equations  (26)  and  (27)  into  equations  (22) 
and  (23)  resulting  in 


a  =  q  -  V  AC.  (t-x)a(x)dx 
1  0  w 


st  t 

+  s5-  j  AC,  (t-x)[d(x)  -  e(x)]dx  -  C  6eV 
0  *  6e  j 

„2C  ft.S,  t 

-  ^21 —  ''“cS~  ;o  AC|_t(t_T)C°l(T)  '  e(x)]dx 


C  ^  C  A 

“V  2U  “  V  °e) 
q  «e 


(28) 


(29) 


Equations  (28)  and  (29)  can  be  solved  for  a  and  q  numerically  or  by  Laplace 
transform  techniques  provided  certain  approximations  for  AC^  and  Ae  are  used 
as  shown  in  the  following  section. 
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2 .  Approximate  Solution  in  Frequency  Domain 

The  time  domain  equations  for  a  and  q  as  given  from  equations  (26)  and 
(27)  are  difficult  to  solve  even  with  modern  high-speed  digital  computers  due 
to  the  presence  of  the  convolution  integrals.  However,  they  are  considerably 
more  tractable  if  they  are  transformed  to  the  frequency  domain  after  making 
the  following  approximations: 

(  biut> 

AC.  =  (C.  )  [1  -a, exp  -  ]  (30) 

w  \  »  '  t  W 


(  b?Ut\ 

acl  (cL  )  [i  -  a2exp  rcTTr]  ^ 

t  at  M  \  til) 

c  c 


where  JC^  |  ,  j  and 

l  awJ  l  atJ 


(|^-j  are  final  steady  state  values  of  AC.  ,  AC.  and 

\d07co  Lw  Lt 


Ac  respectively.  The  constants  a1 ,  bj ,  a 2,  b2>  f ^ ,  g^ ,  and  h-j  are  determined  by 
curve  fits  to  the  results  represented  by  Equations  (16)-(19).  These  approxi¬ 
mations  result  in  the  following  expressions  for  Cz,  and  c: 

Cz  (s)  =  -  sct(s)  ACL  (s)  -  ^  s  (a(s)  -  e(s)  )  ACL  (s)  +  Cz  60(s)  (33) 


Cm  (s)  =  -  s  (5(s)  -  l(s)  )  ACl  (s)  +  Cm  q(s )  +  C  ?  (s)  (34) 

t  q  5e 

e  (s)  =  sa(s)Ae(s)  (35) 


where 


ACl  (s)  =  CL  ]  p/s 

w  l  V 


stbi  f * 


. 1,1  ■  M 


1/S  - 


s+b2  (-T 


(4) 


and 


Ei  is  the  error  integral  whose  series  form  is 

oo  n 

Ei  (s)  =  v  +  £n  s  +  I  -7-t  (39) 

n=l 


v  is  the  Euler  constant  with  numerical  value  0.5772156. 

We  note  that  Cz(s)  and  Cm(s)  can  be  rewritten  in  the  convenient  form 

Cz(s)  =  f (s )  S(s)  +  Cz  6e(s)  (40) 

6e 

Cm(s)  =  g(s)  a(s)  +  Cm  +  Cm  (41) 


Upon  taking  the  Laplace  transforms  of  equations  (21)-(22)  and  substitution  from 
equations  (40)-(41),  the  following  matrix  equation  for  a(a)  and  q(s)  evolves 


from  which  expressions  for  the  transfer  functions 


g(s)  and  q(s)  can  be  found. 
Se(s)  6fi(s) 
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SECTION  V 


APPLICATIONS 

In  order  to  assess  the  effect  of  the  unsteady  terms  on  the  aircraft's 
response  to  control  inputs,  calculations  have  been  made  using  the  Navion 
(Figure  8)  aircraft  as  subject.  This  aircraft  was  chosen  since  its 
characteristics  have  been  studied  extensively  from  flight  research  data  and 
well  documented  by  recent  NASA  technical  reports  [26].  The  geometric  and  mass 
characteristics  of  the  Navion  aircraft  are  listed  in  Table  1  and  the  quasi¬ 
steady  stability  and  control  derivatives  are  listed  in  Table  II.  For  the 
simulation,  a  trimmed  flight  condition  at  5,000  feet  altitude  and  speed  of 
240  ft/sec  was  used.  At  this  time  a  step  input  of  -30°  elevator  deflection 
was  applied  and  held  for  4  seconds.  The  subsequent  responses  of  the  aircraft 
using  both  quasi-steady  and  unsteady  theory  is  shown  in  Figures  9-10.  In 
Figure  9,  the  aircraft's  angle  of  attack  and  pitch  rate  exhibit  a  more  sluggish 
response  to  the  elevator  input  when  unsteady  aerodynamics  are  included.  This  is 
to  be  expected  due  to  the  hysteresis  previously  noted  in  Figure  5.  The  reduction 
in  oscillatory  behavior  in  pitch  rate  is  particularly  noticeable.  From  Figure 
10,  it  can  be  seen  that  the  unsteady  aerodynamics  has  considerably  less  effect  on 
the  aircraft's  forward  speed  and  attitude  response.  These  effects  are  summarized 
in  Table  III. 

Frequency  response  curves  for  the  two  transfer  functions  defined  by  Equation 
(42)  are  shown  in  Figure  11  for  the  Navion  aircraft.  For  comparison  purposes 
the  quasi-steady  form  of  the  transfer  function  is  also  presented.  The  unsteadi¬ 
ness  is  seen  to  have  a  pronounced  effect  on  the  pitch  rate  response  of  the  air¬ 
craft. 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

A  simplified  aerodynamic  model  based  on  the  principle  of  Prandtl's  lifting 
line  theory  and  trailing  vortex  concept  has  been  developed  for  use  in  prelimi¬ 
nary  design  and  stability  and  control  analyses.  The  model  proved  to  be  apolicable 
to  both  two  and  three  dimensional  wing  theory  in  calculating  unsteady  lift  load 
for  arbitrary  angle  of  attack  variation.  Sample  numerical  solutions  to  the 
longitudinal  equations  incorporating  the  unsteady  model  are  presented  for  the 
Navion  aircraft.  From  these  calculations,  an  apparent  increase  in  damping  ef¬ 
fects  in  the  aircraft  response  is  noted. 

For  future  analyses,  it  is  recommended  that  the  control  input  rates  be  in¬ 
cluded  in  the  aircraft  equations  of  motion  to  more  realistically  simulate  and 
assess  the  role  of  the  unsteadiness  in  the  indicial  response.  Further,  a  study 
of  the  aircraft  response  should  be  performed  with  aeroelastic  effects  and  un¬ 
steady  aerodynamics  included. 


TABLE  1 


GEOMETRIC  AND  MASS  CHARACTERISTICS  OF  NAVION 

Aircraft  Mass,  slugs  92.17 

Aircraft  Moment  of  Inertia  (1^),  slug-ft2  —  -  -  -  -  2772.86 

Wing: 

Area,  ft2 . . - . .  .  184.0 

Aspect  Ratio  -------------------  6.04 

Span,  ft -  33.38 

Mean  chord,  ft------------------  5.70 

Horizontal  Tail: 

Area,  ft2  --------  .  -------  43.0 

Aspect  ratio  - -  -  -  - - -  4.0 

Span,  ft  .  12.70 

Mean  chord,  ft------------------  3. 30 

Tail  length  (center  of  gravity  to  tail 

aerodynamic  center),  ft  ------  -  15.12 
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TABLE  3 


STABILITY  CHARACTERISTICS 


Specification 

a 

! 

i 

i 

| 

q 

Quasi -steady 

Unsteady  ' 

i 

Quasi-steady 

Unsteady 

Percent  overshoot 

4.5% 

1% 

85% 

50% 

Time  to  peak 

0.6  sec 

0.9  sec 

0.3  sec 

0.2  sec 

Rise  time 

0.35  sec 

0.45  sec 

0.1  sec 

0.1  sec 

(90%  of  steady 

state  value) 

Settling  time 

1.6  sec 

1.3  sec 

1.3  sec 

0.8  sec 

16 


Equal 

Distance 
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a,  degrees 

Figure  3.  Hysteresis  Effect  Due  to  Unsteadiness 


All  dimensions  are  in  feet  (meters). 

Figure  8.  Sketch  of  Navion  Aircraft 
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oe  aeg  q,  rad/sec 


Figure  9.  Angle  of  Attack  and  Pitch  Rate 
Response  to  Elevator  Step  Input 
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co  (rad/sec) 


Figure  11.  Frequency  Response  of  Navion  Aircraft 


APPENDIX 


VON  KARMAN  -  SEARS  SOLUTION  TO  THE  WAGNER  PROBLEM 

The  problem  Wagner  [4]  solved  was  the  transient  lift  generated  on  a  flat 
plate  airfoil  inclined  to  the  freestream  at  an  angle  of  attack  a  started 
impulsely  from  rest  as  shown  in  Figure  A1 .  The  circulation  around  the  air¬ 
foil  changes  for  unsteadiness  in  either  the  velocity  or  the  angle  of  attack 
of  the  airfoil.  Consequently,  an  equal  amount  of  vorticity  of  opposite  sign 
forms  in  the  wake  and  remains  fixed  relative  to  the  airfoil.  First  we  con¬ 
sider  the  effect  of  an  isolated  free  vortex  of  strength  y(£)d£  located  at  a 
distance  E,  from  the  origin.  The  total  circulation  T  about  the  airfoil  is 
made  up  of  two  parts:  that  due  to  motion  with  no  wake  (quasi-steady  motion) 
with  circulation  Fq  and  the  motion  induced  by  the  wake  which  creates  a  cir¬ 
culation  r1  around  the  airfoil. 

Then 

+1 

r=  rQ  +  =  /  y(x)dx  (Al) 

The  value  of  Tg  (quasi-steady  theory)  is  given  as 

rQ  =  pUc  (A2) 

To  obtain  T-|  we  must  first  determine  y(x),  the  vorticity  distribution 
along  the  plate.  This  is  obtained  easiest  from  methods  of  conformal  mapping 
according  to  Figure  A2.  In  the  z'-  plane,  the  vortex  pair  r*  and  -r*  located 
at  n  and  1/n  respectively  will  create  a  stream  function  given  by 

'f(r>  B)  =  ^  *n  r2  -  Jjj:  £n  r]  (A3) 
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where 


I 


r-|  =  [r2  +  (1/n)2  -2r(l/n)cosB]1//2 


~2  =  [p2  +  n2  -  2rri  cosg]^2 


This  results  in  a  velocity  field  on  the  surface  of  the  cylinder  given  by 


ur  =  0 


II  >1  g  -1 

2tt  C~cosB 


At  the  trailing  edge  of  the  plate  (3  =  0)  we  have 


_  f  U+1 


u3  =  2ir  si  C-l 


Since  we  seek  a  physically  realistic  flow  satisfying  Kutta's  condition 
at  the  trailing  edge  (B  =  0)  we  must  distribute  along  the  plate  a  vorticity 
of  strength  y(x)  which  will  create,  on  the  cylinder,  a  tangential  velocity 


l^l^yand  radial  velocity  zero  so  that  the  sum  of  this  value  and  that 


from  equation  (A5)  results  in  a  surface  tangential  flow  field  satisfying  the 
Kutta  condition  at  the  trailing  edge.  Then  the  resultant  tangential  velocity 
at  the  surface  is  given  as 


UB  =  2tt  l  C-cosB  '  V5-1 


The  corresponding  point  on  the  surface  of  the  flat  plate  has  the  velocity. 


JB  _  r _  r  /g2-i 


sing  2TrsinB  L  ^-cosg  V  1 


Note  that  on  the  surface  of  the  plate,  x  =  cosg.  Also  note  that  on  the  upper 
surface  sing>o  and  on  the  lower  surface  sing<o.  This  results  in 


n,  _  _LH_  /1-x  fe+1 
qu  2tt  C-x  Vl+x  U-l 


Since  y(x)  equals  the  jump  in  tangential  velocity  across  the  plate  we  have 


Y(x)  = 


(All) 


Then  the  circulation  around  the  plate  due  to  the  free  vortex  r'  located 
at  E,  is 


t1 

-1 


y(x)dx  =  — - 


+1 

J 

-1 


1 1-x  dx 
v 1 +x  5-x 


(A12) 


For  a  continuous  distribution  r1  =  y(£)d£,  the  total  circulation  about  the 
plate  due  to  the  wake  is  obtained  by  integrating  over  the  length  of  the  wake 


r-|  =  /w  Ytt)  (J^~  1)  dc  (A13) 

Then  the  total  circulation  about  the  plate  is 

T  -  *Uc  +  /w  y(C)  (  Jp}  -  1)  dC  (A14) 

The  remaining  unknown  is  the  strength  of  the  vorticity  in  the  wake  y(£). 
We  find  this  by  use  of  the  fact  that  the  total  circulation  of  the  whole 
system  (plate  and  wake)  is  zero,  i.e., 

r  +  I  y(c)d£  =  0  (A15) 

1 

or 


ttIIc  + 


£  w 

;  Y(C) 

1 


o 


(A16) 


For  the  impulsively  started  plate  £  =  1+Ut  and  y(0  is  determined  from  the 

w 

integral  equation 
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!  +  1+Ut  xiCl  M  df 

1  ttUc  M  af'  0 

Wagner  solved  this  numerically  for  yCO/wUc. 

The  forces  on  the  plate  can  be  found  from  the  unsteady  form  of 
Bernoulli's  equation 

P  *  -  V2  pq'2  -  ff 

where 

<{>  =  -Uxcosct  -  Uysina  + 

The  forces  in  the  x  and  y  directions  are  then 
X  =  -  ;  ptds 


(A17) 


(A18) 

(A19) 


(A20) 


Y  =  -  ;  pmds 
c 


where  l  and  m  are  direction  cosines  of  the  normal  to  the  surface  and  the 
integration  is  taken  from  the  trailing  edge  over  the  plate  ^nd  back  to  the 
trailing  edge.  Glauert  [5]  has  shown  that  the  integration  reduces  to  the 
following  form: 

X  =  -plirsina 

Y  =  pUTcosa  +  +  p|-jr  ;  q'xdx  (A21) 

Von  KSrnrfn  and  Sears  [9]  compute  the  lift  directly  by  relating  it  to  the 
momentum  of  the  fluid  expressed  as 

L  -  -  (A22) 


where 

I 


y(x)xdx  +  p  j 
1 


(A23) 
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This  results  in  a  sectional  lift  coefficient  given  by  (a  =  1) 
1+Ut  Ac 

ct  ■ 2" ['  *  <  dy 

i  /r- 1 


(A24) 


The  quantity 


1+Ut  tc\  Ac 

Mt)  =  1  +  i  -4P- 

i  ’Uc  /prj 


(A25) 


is  called  the  Wagner  function. 
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Figure  A1 .  Flat  Plate  Started  Impulsively  From  Rest 


z'-plane 
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